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Abstract. We discuss outer and inner crusts of neutron stars in strong magnetic fields. Here, 
we demonstrate the effect of Landau quantization of electrons on the ground state properties 
of matter in outer and inner crusts in magnetars. This effect leads to the enhancement of 
the electron number density in strong magnetic fields with respect to the zero field case. 
For the outer crust, we adopt the magnetic Baym-Pethick-Sutherland model and obtain the 
sequence of nuclei and equation of state (EoS). The properties of nuclei in the inner crust in 
the presence of strong magnetic fields are investigated using the Thomas-Fermi model. The 
coexistence of two phases of nuclear matter - liquid and gas, is assumed in this case. The 
proton number density in the Wigner-Seitz cell is affected in strong magnetic fields through the 
charge neutrality. We perform this calculation using the Skyrme nucleon-nucleon interaction 
with different parameterisations. We find nuclei with larger mass and atomic numbers in the 
inner crust in the presence of strong magnetic fields than those of the zero field case for all 
those parameter sets. Further we investigate torsional shear mode frequencies using the results 
of magnetised neutron star crusts and compare those with observations. 



1. Introduction 

Neutron star crust plays an important role in many observational findings on neutron stars 
[T]. Heat transport and magnetic field evolution in the crust are sensitive to the composition 
of the crust. This, in turn, might influence various transport properties such as electrical and 
thermal conductivities, shear modulus and shear viscosity in the neutron star crust. Figure 
1 demonstrates different regions of a neutron star schematically. The crustal region is just 
beneath the envelope as shown in the figure. Neutron star crust is separated into outer and 
inner crust. The outer crust begins at ~ lO'^g/cm'^. The thickness of the outer crust is a 
few hundred meters. The outer crust is made of nuclei arranged in a body-centered cubic (bcc) 
lattice and embedded in a uniform background of electrons which become relativistic at a density 
~ lO"^ g/cm^. Neutrons and protons are bound in nuclei in this case. The equilibrium nucleus at 
~ W^g/cm^ is ^^Fe. Equilibrium nuclei become more and more neutron rich as density increases. 
The neutron chemical potential exceeds the neutron rest mass at a density ~ 4 x W^^g/cm'^. 
Consequently, neutrons drip out of nuclei and it is called the neutron drip point. This is the 
end of the outer crust and the beginning of the inner crust. The inner crust is composed of 
neutron rich nuclei arranged in a bcc lattice and also in coexistence with dripped neutrons and 
a uniform background of electrons. The thickness of the inner crust is ~ 1 km. More and more 
neutrons come out of neutron rich nuclei with increasing density. Finally, nuclei of the inner 



crust dissolve into a uniform nuclear matter at the crust-core boundary. 
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Figure 1. Schematic structure of a neutron star interior. 



Extreme physical condition such as very strong surface magnetic fields ~ 10 G were found 
to exist in a class of neutron stars called magnetars [21 [3]. Even stronger interior field could 
exist in magnetars. Magnetars are classified into two groups - Anomalous x-ray pulsars (AXPs) 
and soft gamma-ray repeaters (SGRs) [H [5]. SGRs are found to emit sporadic bursts of soft 
gamma rays. It was observed that SGRs occasionally emit stronger gamma rays known as giant 
flares having luminosities as high as 10^^ ergs s~^ O 13 13 [9] . Giant flairs are considered to be 
the result of evolving magnetic fleld and its stress on the magnetar crust. In the decaying tail of 
giant flairs, quasi-periodic oscillations (QPOs) were observed over long duration [71|8l|9]. QPOs 
were identified as torsional shear modes of magnetar crust [3]. The study of QPOs might shed 
light on the ground state properties of magnetar crust. 

Recent studies showed that high magnetic fields associated with magnetars would influence 
the composition and EoS of neutron star crusts which, in turn, might have impact on the QPO 
frequencies. In this paper, we study and discuss this problem. The paper is organised in the 
following way. We discuss the composition and EoS of neutron star crusts with or without 
magnetic field in section 2. Torsional shear modes are discussed in section 3. We summarise in 
section 4. 

2. Composition and EoS of Neutron Star Crust 

In a seminal paper, Baym, Pethick and Sutherland (BPS) studied the equilibrium composition 
and EoS of the ground state matter of the outer crust without magnetic field [10]. Later this 
BPS model was extended to include the effects of strong magnetic fields [llj. High magnetic 
fields of magnetars might infiuence the composition and EoS through Landau quantisation. 
Recently we have revisited the magnetic BPS model. Furthermore, we include the finite size 
effect in the lattice energy and adopted recent experimental and theoretical mass tables jl2j. We 
adopt the Wigner-Seitz approximation in this calculation. The equilibrium sequence of nuclei 
is determined by minimising the Gibbs free energy per baryon (g) with respect to mass number 
(A) and atomic number (Z) at a constant pressure (P) |10] . 
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where Ub is the baryon number density and tin is the number density of nuclei; those two are 
connected by Ub = Atln. The energy of the nucleus (including rest mass energy of nucleons) is 



Wn = rriniA - Z) + rUpZ - bA , (2) 

where b is the binding energy per nucleon. We obtain experimental nuclear masses from the 
atomic mass table of Ref.([13j). We exploit the theoretical extrapolation for the rest of nuclei 
[H] . The lattice energy of the cell Wl is given below [15] , 



The cell and nucleus radii are rc and r^, respectively. The first term in Eq.(l3]) is the lattice 
energy for point nuclei and the second term corresponds to the finite size effect [J^. 

Now we focus on the impact of strong magnetic fields on the Gibbs free energy per baryon. In 
presence of a magnetic field, the motion of charged particles perpendicular to the magnetic field 
is quantised into Landau levels. This is known as Landau Quantisation [16j . In this calculation, 
we are interested only in Landau quantisation of electrons and its infiuence on ground state 
properties of neutron star crusts. The electron energy density (eg), pressure (Pe) and number 
density (ng) due to Landau quantisation of electrons are given below [12j . 
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where B is the magnetic field, Umax is the maximum Landau quantum number, pp^ ^ = 
\/ fj-'i — {mi + 2i>eB), /ig is the electron chemical potential, and the spin degeneracy 51,^ = 1 
for = and = 2 for other values of Landau quantum numbers. We define B^ = B/Bc, 
where the critical field for electrons is Be = 4.414 x 10^'^ G. In the following paragraphs, we use 
magnetic field strengths in terms of B^. 

Electron number density is greatly influenced due to Landau quantization in strongly 
quantising magnetic fields. For B^ < 10^, electrons populate large number of Landau levels. 
Consequently, there is no change in electron number density compared with the zero field 
situation. However, electron number density strongly enhanced when the field is -B > 10^^ 
G. This is attributed to the fact that electrons populate only the zeroth Landau level [I^. This 
might impact the composition and EoS of the outer crust. 

We calculate the sequence of nuclei in the outer crust for i? = and B^ = 10^. Results 
are shown in Table 1 and Table 2. A careful examination of results in both tables shows that 
nucleus such as ^^Ni which is present for = 0, is absent in the case of B^ = 10^ Similarly, 
two new nuclei ^^Sr and ^^^Pd are found to appear in Table 2. The finite size correction to the 
lattice energy modifies the sequence of nuclei for 5* = 10^. It is to be noted that our results are 
different from those of Ref.[TT] because we have used the recent experimental and theoretical 
nuclear mass tables in our calculation. One important result of this calculation is the shift of the 
neutron drip point to higher density in presence of strong magnetic fields. This is evident when 
we compare pmax which is the maximum density for the existence of an equilibrium nucleus, in 
the last row of Table 1 and Table 2. 



Table 1. Sequence of nuclei in the outer crust in absence of magnetic fields. Elements with 
their charges (Z) and neutron numbers (A^), maximum mass density (pmax), electron number 
density (ng), electron chemical potential (/j,) and Gibbs free energy per nucleon (g) for 5 = 
are recorded here. 



element 


Z 


N 


Pmax 






n 


e 


Pe 


g 








(g/cm^) 




(cm 




(MeV) 


(MeV) 




26 


30 


8.01 X 10'* 


2 


24 


X 




0.95 


930.60 


62Ni 


28 


34 


2.71 X 10^ 


7 


38 


X 


1031 


2.61 


931.32 


64Ni 


28 


34 


1.33 X 10^ 


3 


51 


X 


10^2 


4.34 


932.04 


66Ni 


28 


34 


1.50 X 10^ 


3 


82 


X 


1032 


4.46 


932.09 


86Kr 


36 


50 


3.10 X 10^ 


7 


80 


X 


1032 


5.64 


932.56 


84Se 


34 


50 


1.06 X 10^" 


2 


58 


X 


1033 


8.39 


933.62 


82 Gc 


32 


50 


2.79 X 10^" 


6 


54 


X 


1033 


11.43 


934.75 


80Zn 


30 


50 


6.11 X 10^° 


1 


37 


X 


103'1 


14.63 


935.90 




28 


50 


9.25 X 10^° 


1 


99 


X 


1034 


16.56 


936.57 


126Ru 


44 


82 


1.29 X 10^^ 


2 


69 


X 


1034 


18.30 


937.12 


124Mo 


42 


82 


1.86 X 10^^ 


3 


78 


X 


1034 


20.50 


937.83 


1222J. 


40 


82 


2.64 X 10^^ 


5 


18 


X 


1034 


22.76 


938.53 


i20Sr 


38 


82 


3.77 X 10" 


7 


13 


X 


1034 


25.33 


939.31 


iiSKr 


36 


82 


4.34 X 10" 


7 


91 


X 


1034 


26.22 


939.57 



Table 2. Same as Table 1 but for 5* = 103. 



element 


Z 


N 


Pmax 




n 


e 


Pe 


g 








(g/cm-^) 




(cm" 




(MeV) 


(MeV) 




26 


30 


4.31 X 


10^ 


1 


21 X 


1033 


0.87 


930.46 


62Ni 


28 


34 


1.83 X 


lOio 


4 


99 X 


1033 


2.94 


931.32 


64Ni 


28 


34 


2.33 X 


lO^o 


6 


14 X 


1033 


3.60 


931.59 


88Sr 


38 


50 


2.59 X 


10^0 


6 


73 X 


1034 


3.94 


931.72 


86 Kr 


36 


50 


4.33 X 


lOio 


1 


09 X 


1034 


6.35 


932.69 


84Se 


34 


50 


6.33 X 


lOio 


1 


54 X 


1034 


8.97 


933.72 


82Ge 


32 


50 


8.67 X 


IQio 


2 


04 X 


1034 


11.83 


934.82 


80Zn 


30 


50 


1.13 X 


10" 


2 


55 X 


1034 


14.82 


935.93 


128pd 


46 


82 


1.29 X 


10" 


2 


78 X 


1034 


16.15 


936.38 


126Ru 


44 


82 


1.50 X 


10" 


3 


15 X 


1034 


18.32 


937.13 


124Mo 


42 


82 


1.72 X 


10" 


3 


50 X 


1034 


20.35 


937.82 


i22Zr 


40 


82 


1.96 X 


10" 


3 


86 X 


1034 


22.45 


938.50 


120gj. 


38 


82 


4.34 X 


10" 


8 


22 X 


1034 


25.44 


939.32 




3G 


82 


4.92 X 


10" 


8 


98 X 


103 i 


20.27 


939.57 



Now we discuss the inner crust. Nuclei of the inner crust are not only immersed in a uniform 
background of electrons but also in coexistence with a neutron gas. Besides charge neutrality, 
the matter of the inner crust is in /3-equilibrium. The ground state properties of the inner crust 
in zero magnetic field were studied by several groups |171 [181 [IHl EQl Elj • The first calculation of 
the inner crust in presence of strong magnetic fields was investigated in the Thomas-Fermi (TF) 
model at zero temperature [2^. Nuclei of the inner crust are also arranged in a bcc lattice. The 
Wigner-Seitz approximation is used in the calculation of the inner crust in presence of strong 
magnetic fields. Further, each lattice volume is replaced by a spherical cell with one nucleus at 
the center. Each cell is charge neutral. There is no Coulomb interaction between cells. In this 
case, the denser phase in a nucleus is in coexistence with the low density phase of neutron gas. 
This has to be dealt with in a thermodynamically consistent way. The spherical cell in which 
protons and neutrons reside, does not define a nucleus. We exploit the subtraction procedure 
of Bonche, Levit and Vautherin (BLV) |23] to remove the gas part from the nucleus plus gas in 
the cell. In the BLV scheme, the density profiles of the gas part as well as nucleus plus gas are 
derived self-consistently in the TF formalism. Finally, the difference of two solutions in BLV 
scheme determines the nucleus \23\ 123] . Here electrons are Landau quantised in strong magnetic 
fields, but protons in the cell are indirectly affected through the charge neutrality. 

Our starting point in this calculation of inner crust is the thermodynamical potential (Oat) 
of a nucleus. This is obtained as the difference of the thermodynamical potential in nucleus plus 
gas {^Ing) and that of the gas phase (J^g)- This is given by 



Otv = ^NG — ^G , (5) 

where the thermodynamic potential is written in terms of free energy (F), chemical potential 
(Hi) and number density (rij) of i-th species, 

r2 = F - ^ fiiUi . (6) 

i=n,p 

The free energy is a function of average baryon density (rib) and proton fraction (1^) and given 
below, 

F{nb, Yp) = J[eN + ec + ee]dr . (7) 

The free energy contains nuclear energy density (en), Coulomb energy density (Sc) and electron 
energy density. We calculate the nuclear energy density using Skyrme (SkM) nucleon-nucleon 
interaction [22 | 125 1 l26 l I27j. We minimise the thermodynamical potential with respect to baryon 
density in each phase and obtain neutron and proton density profiles of two phases solving 
coupled equations [22l[28]. The chemical potential in presence of magnetic fields is [22j 



Me 



■ o o 1 1/2 

pk^^+ml + 2eBu - < y-(r) > , (8) 



where < V'^{r) > denotes the average single particle Coulomb potential. As soon as we know 
density profiles in two phases, we calculate mass and atomic numbers of a nucleus using the 
subtraction procedure of BLV [22] . 

Now we present our results of the inner crust. As electrons are Landau quantised in strong 
magnetic fields, we show its impact on electron number density. Figure 2 shows electron number 
density as a function of mass density. For B^^ = 10^ or smaller field strength, large numbers 
of Landau levels are populated by electrons. As a result, the electron number density does 
not differ from the zero field results. For stronger magnetic fields such as = 10^, electrons 
populate the zeroth Landau level or a few levels over a certain mass density regime. In this 
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Figure 2. Electron number density is plotted with mass density for different magnetic field 
strengths. 

case, it is evident from Fig. [2] that electron number density is significantly enhanced compared 
with the zero field case. However, the situation is different at higher densities. We do not see 
any enhancement of electron fraction there. Proton fraction is modified in presence of strong 
magnetic fields through charge neutrality. Neutrons are affected by magnetic fields through 
/3-equilibrium. The outcome of this is fewer neutrons drip out of a nucleus in presence of strong 
magnetic fields |22j . 

Next we obtain the equilibrium nucleus at each density point by minimising the free energy 
per particle of the cell [22]. In the left panel of Fig. [3l we show the sequence of equilibrium 
nuclei and their corresponding free energies per particle as a function of mass density with and 
without magnetic field. Furthermore, this calculation is done without subtracting the gas part. 
It is noted that the free energy per nucleon is reduced in strong magnetic fields of B^, = 10^. This 
is attributed to the population of electrons in the zeroth Landau level. However, free energy per 
nucleon for B = and B^, = 10^ come closer when large number of Landau levels are populated 
at higher densities. For zero magnetic field case, we find excellent qualitative agreement between 
our results and those of Negele and Vautherin [20j . In the right panel of Fig. [3l mass and atomic 
numbers of equilibrium nuclei are shown after applying the subtraction method of BLV. We find 
that mass and atomic numbers increase in presence of a strong magnetic field compared with 
the results of zero magnetic field. 

We repeat the calculation of inner crust with SLy4 [29j and Sk272 [30j nucleon-nucleon 
interactions. Results of these calculations are shown in Fig. [H We find qualitatively similar 
results for mass and atomic numbers as a function of average baryon number density as found 
in case of SkM interaction. For SLy4 interaction, we find mass and atomic numbers in the 
magnetic field jump when electrons in the zeroth Landau level move to the next level around 
baryon number density 0.05 fm~^. Further, we note that the SLy4 parameter set results in 
higher mass and atomic numbers because of stiffer density dependence of the symmetry energy 
at sub-saturation densities than that of other parameter sets considered here. 

3. Torsional Shear Modes 

Next we are interested in the role of compositions and EoS of magnetised crusts that we have 
described in section 2, on torsional shear modes of magnetars. It is argued that torsional shear 




Figure 3. Free energy per nucleon of the cell is plotted with mass density for zero magnetic 
field, B = 10^ and 5* = 10^ (left panel). Equilibrium nuclei are denoted with solid symbols in 
both panels; mass and atomic numbers in nuclei after the subtraction of the gas phase is shown 
as a function of average baryon number density (left panel). 
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Figure 4. Mass and atomic numbers in nuclei after the subtraction of the gas phase is shown 
as a function of average baryon number density for SLy4 interaction (left panel) as well as Sk272 
(right panel) with B = and B^, = 10^. 



mode frequencies depends on the shear modulus of neutron star crusts [HI [9l [3T]. Further the 
shear modulus is sensitive to compositions and EoS of neutron star crusts. This implies that 
nuclear physics plays an important role in determining torsional shear frequencies j31| . Here 
we describe the calculation of the shear modulus and shear velocity of the magnetised neutron 
star crusts. We use the following the expression of the shear modulus(/u) at zero temperature 

MM, 

^ = 0.1194^^ , (9) 

a 

where a = 3/(47rnj), Z is the atomic number of a nucleus and rii is the ion density. This form of 
the shear modulus was obtained by assuming a bcc lattice and performing directional averages 



|34j . We exploit inputs rij and a from the magnetised crusts of section 2. We immediately 
calculate the shear speed as Vs = ^/{p/p). In Figure 5, the shear speed is shown as a function of 
mass density with and without magnetic field. For B^^ = 10^, the shear speed does not differ from 
that of -B = 0. In presence of strong fields such as B^: = 10^, the shear speed enhanced compared 
with the zero field results. We find spikes in each curve. It happens when an equilibrium nucleus 
jumps to a new one in the outer crust. 
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Figure 5. Shear velocity is shown as a function of mass density for different magnetic field 
strengths. 

Torsional shear modes were studied extensively using non-magnetic crusts in Newtonian 
gravity [3l [351 EH EZ] and general relativity [381 ESI HOI HH US]. Here we are interested in 
the effects of magnetised crusts on torsional shear mode frequencies. We adopt the model of 
Refs. [38ll42| to calculate torsional shear frequencies in presence of a dipole magnetic field. In this 
calculation, we consider the magnetised crust decoupled from the core. Further, the magnetised 
star is considered to be the spherically symmetric because the deformation due to the magnetic 
field is neglected |42j. The perturbed equations describing torsional shear modes are obtained 
by linearising the equations that govern equilibrium configurations. Equilibrium stellar models 
are determined using the metric 

ds'^ = -e'^^dt'^ + e^^dr^ + (dO^ + sin'^edcj)'^) . (10) 

We consider axial-type perturbation in the four velocity {u^^) and the only non-vanishing 
perturbed quantity is the 4> component of the perturbed four velocity w'^ [38] 

du't' = e-'^dty{t, r)^dePi{cose) , (11) 

where dt and de correspond to partial derivatives with respect to time and 9, respectively, 
Pi{cos9) is the Legendre polynomial of order I and y{t,r) is the angular displacement of the 
matter. Assuming a harmonic time dependence of y{t, r), one arrives at the eigenvalue equation 
which is a second order differential equation [38] . Finally, we estimate eigenfrequencies by solving 




Figure 6. Fundamental {n = 0, i = 2) shear mode frequency as a function of normalised 
magnetic field for a neutron star of 1.4 Mq (left panel); Torsional shear mode frequency (n = 0) 
as a function of £ values for different magnetic field strengths and a neutron star of 1.4 Mq 
(right panel). 

two first order differential equations Eq.(69) and (70) of Sotani et al. with appropriate boundary 
conditions [38] . 

Now we discuss the results of our calculation using the magnetised crusts calculated with SkM 
nucleon-nucleon interaction. In the left panel of Figure 6, fundamental shear mode frequency 
(n = 0, £ = 2) is plotted with magnetic field for a neutron star of 1.4 Mq. The fundamental 
frequency increases with increasing magnetic field. This plot implies if we know the mass of the 
star and the field strength accurately, we can determine the frequency and compare it with the 
observed frequency. In the right panel of Figure 6, we show torsional shear mode frequency for 
n = as a function of i values. As i values increase, the frequency increases. It is also evident 
from this figure that the torsional frequency increases with increasing magnetic field strengths. 
It is also noted that frequencies of torsional shear modes decrease with increasing masses of 
magnetars |43j . 

Finally, we calculate fundamental torsional shear mode frequencies as well as various overtones 
for SGR 1806-20 and SGR 1900+14 and compare those with observed frequencies of QPOs. 
We show this comparison in Table 3 and Table 4. In case of SGR 1806-20, higher torsional 
frequencies i.e. (93, 150, 626 and 1838 Hz) are nicely explained by our theoretical calculation 
with magnetised crusts when the mag netic field is B = 8 x 10^"^ G and mass is 1.4 Mq [43j . 
However, the agreement of three lower frequencies (18, 26 and 29 Hz) with our calculation is 
not good. These frequencies could possibly be explained by the Alfven modes [44J. On the 
other hand, all four observed frequencies (28, 55, 82, 154 Hz) in SGR 1900+14 are in excellent 
agreement with our calculation when the mass is 1.2 Mq and magnetic field B = 4 x 10^^ G 

m- 

4. Summary 

We have developed a model of magnetised neutron star crusts. In this context, we have 
investigated the effects of Landau quantisation of electrons on the ground state properties 
of outer and inner crusts. It is observed that composition and EoS of neutron star crusts 



Table 3. Calculated torsional shear mode frequencies are compared with observed frequencies 
of SGR 1806-20. 



SGR 1806-20 


Observed 


Calculated 






Frequency 


Frequency 


n 


I 


(Hz) 


(Hz) 






18 


15 





2 


26 


24 





3 


29 


32 





4 


93 


93 





12 


150 


151 





20 


626 


626 


1 


29 


1838 


1834 


4 


2 



Table 4. Same as Table 3 but for SGR 1900-M4. 



SGR 1900+14 


Observed 


Calculated 






Frequency 


Frequency 


n 


I 


(Hz) 


(Hz) 






28 


28 





4 


54 


55 





8 


84 


82 





12 


155 


154 





23 



are appreciably modified in strong magnetic fields. We have calculated the shear modulus 
and shear velocity in connection with torsional shear modes. It is found that torsional shear 
frequencies calculated using our model of magnetised crusts are in good agreement with the 
observed frequencies of QPOs. 
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